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Quantum light-matter interfaces (QLMIs) con¬ 
necting stationary qubits to photons will enable 
optical networks for quantum communications, 
precise global time keeping, photon switching, 
and studies of fundamental physics. Rare-earth- 
ion (REI) doped crystals are state-of-the-art ma¬ 
terials for optical quantum memories and quan¬ 
tum transducers between optical photons, mi¬ 
crowave photons and spin waves. Here we demon¬ 
strate coupling of an ensemble of neodymium 
REIs to photonic nano-cavities fabricated in the 
yttrium orthosilicate host crystal. Cavity quan¬ 
tum electrodynamics effects including Purcell 
enhancement (F=42) and dipole-induced trans¬ 
parency are observed on the highly coherent 
4 -^ 9/2 — 4 ^ 3/2 optical transition. Fluctuations in the 
cavity transmission due to statistical fine struc¬ 
ture of the atomic density are measured, indi¬ 
cating operation at the quantum level. Coher¬ 
ent optical control of cavity-coupled REIs is per¬ 
formed via photon echoes. Long optical coher¬ 
ence times (X 2 ^100 /is) and small inhomogeneous 
broadening are measured for the cavity-coupled 
REIs, thus demonstrating their potential for on- 
chip scalable QLMIs. 

QLMIs are quantum devices composed of light emitters 
with quantum states that can be controlled via optical 
fields and entangled to photons mm- They enable dis¬ 
tribution of quantum entanglement over long distances 
in optical quantum networks for quantum communica¬ 
tions pQ. Quantum networks of atomic clocks have also 
been proposed for precise global time-keeping and stud¬ 
ies of fundamental physics |2|. Realizing QLMIs requires 
control of light and matter at the single atom and sin¬ 
gle photon level, which enable optoelectronic devices like 
optical modulators and nonlinear optical devices operat¬ 
ing at the most fundamental level [3]. QLMIs are also 
expected to play a leading role in realizing optical to 
microwave quantum transducers for interconnecting fu¬ 
ture superconducting quantum machines via optical fi¬ 
bres [4j[5]. 

Scalable and robust QLMIs require emitters to have 
long spin coherence times and coherent optical transi¬ 
tions. For integrated optical quantum networks, these 
emitters need to be coupled to on-chip optical resonators 
that capture the photons in a single mode and further 
couple them into optical fibres or waveguides. The solid- 
state emitters most investigated so far for on-chip QLMIs 


are semiconductor quantum dots (QDs) [B and nitrogen 
vacancy centers in diamond (NVs) [7j. To date, com¬ 
plete quantum control of single QD and NV spins, spin- 
photon entanglement, and entanglement of remote NVs 
via photons have been realized [8HTQ]. Both QDs m 
and NVs pQ have been coupled to optical nano-cavities. 
However, the challenge in growing optically identical QDs 
limits their prospects for a scalable architecture 0. NVs 
embedded in nanostructures have long electronic spin co¬ 
herence times m , but suffer from optical spectral insta¬ 
bilities such as blinking and spectral diffusion [14]. These 
spectral instabilities have so far impeded the coherent 
coupling between optical fields and NV centers in nano¬ 
resonators that are essential for further developments of 
QLMIs. 

Rare-earth ions (REIs) embedded in host crystals at 
cryogenic temperatures exhibit highly coherent quantum 
states in the 4f orbital m- The Zeeman or hyper- 
fine states of REIs can have coherence times as long 
as six hours [16], the longest ever demonstrated in a 
solid. These states are connected via optical transitions 
with the narrowest linewidth in the solid state (sub-kHz) 
and small inhomogeneous broadening (MHz to GHz) [17] . 
This outstanding optical and spin coherence makes REI- 
doped crystals the state-of-the-art material for macro¬ 
scopic solid-state optical quantum memories [181 [19] . In¬ 
tegrated REI-doped waveguide quantum memories have 
also been developed 123 EI]. Detection and control of sin¬ 
gle REI spins has been recently demonstrated in bulk ma¬ 
terial, but not using the transitions employed in optical 
quantum memories at cryogenic temperatures [22, [23] . 
Coupling the highly coherent optical transitions of REIs 
to nano-cavities will enable on-chip QLMIs where REI 
ensembles act as quantum memories and single REIs act 
as qubits [24]. 

Here we demonstrate high-cooperativity coupling of a 
neodymium (Nd 3+ ) ensemble to photonic nano-cavities 
fabricated directly in the yttrium orthosilicate (YSO) 
host crystal, and show coherent optical control of REIs 
coupled to nanophotonic cavities. These results are en¬ 
abled by the long coherence time and small inhomoge¬ 
neous broadening of cavity-coupled REIs, which are es¬ 
sential properties that may lead to nanophotonic QLMIs 
with better prospects for scalability than those based on 
NVs and QDs. 
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FIG. 1. Photonic crystal nano-beam resonator fabricated in Nd:YSO. a. Scanning electron microscope image of the 
device. The red inset is a zoomed-in view of the 45° angle-cut coupler that allows vertical coupling of light from a microscope 
objective. The blue inset shows the grooves forming the photonic crystal, b. Schematics of the nano-beam resonator with 
simulated field profiles of the fundamental TE resonance mode. The TE polarization aligns with the Di axis of the YSO crystal. 
A magnetic field of 500 mT is applied in the D 1 -D 2 plane at an angle of a=135° with respect to Di axis. c. Broadband cavity 
transmission spectrum showing the cavity resonance with quality factor Q=4,400. 


RESULTS 

Design, fabrication and characterization of nano¬ 
cavities. The nano-cavities, one of which is shown 
in Fig. [lja), were fabricated in neodymium-doped yt¬ 
trium orthosilicate (Nd 3+ :YSO) using focused ion beam 
milling. For this study, we used devices fabricated 
in two types of samples with Nd doping of 0.2% and 
0.003% (Scientific Materials Inc.). The photonic crystal 
cavity is made of grooves milled in a triangular nano¬ 
beam [25] (Fig. [ljb)) (see Methods). Finite-difference 
time-domain (FDTD) simulations [26] indicate a TE 
mode with quality factor exceeding lxlO 5 , mode vol¬ 
ume Vmode = 1.65(A/nyso) 3 = 0.2/im 3 , and mode pro¬ 
file shown in Fig. [ljb). Here V is defined as U mo de = 
f v e(r)|E(r)| 2 d 3 r/max(e(r)|E(r)| 2 ), where E(r) is the elec¬ 
tric field and e(r) is the electric permittivity at position 
r. Two 45° angled cuts at both ends of the nano-beam 
(i.e. couplers) allow for coupling light from a direction 
normal to the chip (i.e. b in Fig. Il[b)) using a confo- 
cal microscope setup (See Methods). A broadband light 
source was coupled into the resonator from one end and 
the transmitted light was collected from the other coupler 
with typical efficiencies ranging from 20% to 50%. The 
transmitted spectrum shows a resonance with quality fac¬ 
tor Q=4,400 (Fig.[l|c)) in the device used for the follow¬ 
ing measurements. Arrays of devices were reproducibly 
fabricated with similar performance (Supplementary SI). 

Measurement of the coupling rate between REIs 
and the nano-cavity. The coupling of Nd 3+ ions to 
the nano-cavity was observed through enhancement in 


photoluminescence (PL) and emission rates. With the 
0.2% device cooled at 3.5 K, an 810 nm laser coupled 
into the cavity excited PL in the 4 /g/2 — 4 P3/2 tran¬ 
sition that was then collected from the output coupler 
(Fig-H a)). The PL spectrum shows two lines at 883.05 
nm and 884.06 nm, corresponding to two inequivalent 
sites (Yi and Y 2 ) of Nd 3+ in YSO. An important obser¬ 
vation is that the inhomogeneous linewidth of the ions in 
the cavity is the same as in the bulk material, for both the 
0.2% (A in hom=16.0 GHz) and 0.003% (5.9 GHz) devices 
(Supplementary S2). A small inhomogeneous linewidth 
(on the order of ^ 10 GHz) is important for scaling to 
networks of multiple QLMIs (Supplementary S3). The 
cavity resonance was tuned across the Nd 3+ PL line us¬ 
ing a gas condensation technique [F. The spectrograms 
in Fig. [2|b-d) show enhancement of the Yi line when 
the cavity is resonant with it, a signature of coupling. 
The Y 2 line exhibits negligible enhancement because its 
dipole moment is not aligned with the TE cavity polar¬ 
ization (Di axis of the YSO crystal (Fig. 0b)). The 
spontaneous emission rate enhancement was character¬ 
ized via lifetime measurements. A pulsed laser at 810 
nm excited fluorescence of the Yi line, which was filtered 
using a monochromator and detected with a single pho¬ 
ton counter (Fig. [2je)). From single exponential fits, we 
calculated a reduction in lifetime from 254 jus when the 
cavity is detuned by AA =0.3 nm, to 87 jus on resonance. 
Taking into account the branching ratio of the 883 nm 
transition (j8 ~ 4.5%, see Methods), the reduction in life¬ 
times corresponds to an ensemble averaged Purcell fac¬ 
tor [2] F=42, which agrees well with the estimations that 
assume a uniform spatial distribution of Nd 3+ ions in the 
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FIG. 2. Coupling of Nd 3+ ions to the YSO cavity mode showing enhanced photoluminescence (PL) and reduced 
lifetimes, a. Schematic of energy levels for Nd 3+ in YSO. Optical excitation at 810 nm results in PL at several wavelengths 
with only the 883 nm transition enhanced by the cavity, b. Spectrogram showing the Nd 3+ PL while the cavity is tuned 
across resonance using gas condensation. The dashed line is a guide to the eye indicating the central wavelength of the cavity 
resonance. The cavity resonance is not visible because there is no background luminescence to populate the cavity mode. c. 

d. PL spectra in the uncoupled (c) and coupled (d) cases. The cavity resonance was drawn to indicate the cavity location. 

e. Lifetime measurements for coupled (r c =87 /xs, AA = 0) and uncoupled (t°=254 /xs, AA=0.3 nm) cases, f. Change in 
lifetime as a function of the cavity detuning, which fits well with the calculation (red curve) using quality factor Q=4,400, 4.5% 
branching ratio and field intensity averaged over the mode volume. 


resonator (Supplementary S4). A single ion positioned 
at the maximum cavity field would experience a Purcell 
factor ~200. A similar result was obtained in a 0.003% 
cavity with lower quality factor (Supplementary S5). 

Measurement of the optical coherence time for 
cavity-coupled REIs. Coherent and stable optical 
transitions are essential for QLMIs. We characterized the 
optical coherence time T 2 of the 883 nm transition using 
two-pulse (tt/2 - 7r) photon echo techniques (Fig. §a)), 
with an applied magnetic field of B=0.5 T (see Methods). 
The laser pulses were coupled in and the echoes were col¬ 
lected via the couplers while the 0.2% and 0.003% cav¬ 
ities were on resonance with the Nd transition. Since 
only a small sub-ensemble (<100 ions) in the cavity was 
excited, the weak echo signal required detection using 
single photon counters. A typical echo from the 0.2% 
cavity is shown in Fig. [3^d) . The echo decays as a func¬ 
tion of the (7r/2 - 7 r) time delay r are plotted in Fig.J^b) 
together with photon echoes from bulk substrates. For 
the 0.2% sample, T 2 =4r c =2.8d=0.4 /xs was measured for 
the cavity, which shows a good agreement with the bulk 
value of T|=3.2±0.4 /xs. For the 0.003% sample, the 
echo exhibited two exponential decays. The slower de¬ 
cays give T 2 =100±5/xs (bulk) and T 2 =94±5 /xs (cavity), 
which match with values reported in |3]. The fast de¬ 
cays are likely due to the superhyperfine interactions be¬ 


tween Nd 3+ and its neighbouring yttrium ions, which 
commonly manifests as modulated echoes decaying faster 
than T 2 [31. No oscillations was observed in Fig. [3ja) be¬ 
cause of the fast modulation frequency (^1 MHz) due to 
our strong magnetic field. By lowering the B field, we 
observed oscillations of echoes for the initial 10 /xs de¬ 
lay. Additionally, we varied the excitation powers and 
did not see changes in the T 2 values, which indicates our 
measurement was not significantly affected by instanta¬ 
neous spectral diffusion. In sum, the good agreement on 
T 2 between the cavity and bulk confirms that the optical 
coherence property of Nd 3+ ions is not affected by the 
nano-fabrication. For higher Purcell factors, the T 2 in 
cavities should decrease due to the T 2 < 2Ti limit and 
would become smaller than the bulk value. This regime 
is not reached here because the Purcell enhanced 2Ti is 
not smaller than T 2 . 

The observation of photon echoes demonstrates coher¬ 
ent optical control of the quantum state of cavity-coupled 
ions. This control was further extended by varying the 
7r pulse duration and observing Rabi oscillations in the 
echo intensities as shown in Fig. j3^c). A Rabi frequency 
of ^6 MHz is estimated. The same oscillation was not 
observed in the bulk. For the coupled laser power, the op¬ 
timal 7 r pulse duration is 0.4 /xs. The oscillations are not 
visible for pulse duration less than 0.3 /xs, because of the 
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FIG. 3. Photon echo measurements showing long optical coherence time and enhanced echo intensity from an 
ensemble of Nd 3+ ions in the cavity, a. Two-pulse photon echo sequence ( tt /2 — 7r) used to measure Th. b. Two-pulse 
photon echo decays measured in both the cavity (red) and the bulk (black) samples with two different doping concentrations. 
The inset plots the echo decays measured with a 0.2% doped sample, c. Oscillation of echo intensity with increasing width 
of the 7r rephasing pulse. The periodic signal reveals the ensemble averaged Rabi frequency of the coupled ions. The ideal i r 
pulse duration for the input power was 0.4 g s. d. Enhanced photon echo intensity (by ~12 fold) when the cavity is coupled, 
compared to the uncoupled case (cavity detuned by AA=15 nm so that the transition is outside the photonic bandgap). 


limited rise/fall times (^ 200 ns) of the pulse-generating 
setup (Methods). A ^12 fold increase in the echo inten¬ 
sity is observed in the cavity-coupled case compared to 
the uncoupled case (i.e. detuning AA=15 nm) as shown 
in Fig. §d). This enhancement can be attributed to a 
combination of several effects: the higher atomic absorp¬ 
tion rate through the Purcell effect [35], stronger intra¬ 
cavity field intensity, and high echo collection efficiency 
as the ions emit dominantly into the cavity mode. The 
spectral diffusion of the coupled ions using three-pulse 
photon echoes was also investigated. The homogeneous 
linewidths were broadened at rates of 6.1 kHz //s -1 for 
the 0.2% doped cavity and 380 Hz /is -1 for the 0.003% 
cavity. These slow spectral diffusions permit repeated op¬ 
tical addressing of the ions for 10 s of gs (Supplementary 
S 6 ). 


provided by nano-resonators is that efficient atom-photon 
interaction can be achieved in a small volume with only 
a handful of ions. This is readily observable in our sys¬ 
tem, where the Nd 3+ ions coherently interact with the 
intra-cavity field and control its transmission via dipole- 
induced transparency [29] . With the cavity tuned to 883 
nm, the cavity transmission was probed using broadband 
light and a dip was observed at resonance (Fig. EM. 
The depth of the dip depends on the collective coupling 
cooperativity 77 = 47V g 2 /(ftTh), where g is the ensemble 
averaged coupling strength, k is the cavity full linewidth, 
Th is the Nd 3+ homogeneous linewidth, and N is the 
number of ions per Th- Considering an empirical collec¬ 
tive dipole-cavity coupling model, the normalized cavity 
transmission in the presence of unsaturated resonant ions 
is, 


Measurement of dipole-induced transparency and 
statistical fine structure. QLMIs require efficient in¬ 
teractions between atoms and photons, which is why 
quantum memories use long atomic clouds or doped crys¬ 
tals to achieve large optical depth. One key advantage 


T = 


K 

iA + k + 47V<7 2 /7 h 


2 


(i) 


which simplifies to T = (1 + 77 ) 2 for zero detuning [30] . 
The cavity transmission can be controlled by varying the 
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FIG. 4. Control of cavity transmission and observation of statistical fine structure (SFS) of coherently driven 
Nd 3+ ions in the cavity, a. Broadband transmission spectra as the cavity is tuned to the 883 nm Nd transition. A dip is 
observed when the two are on resonance. The negligible dip at far detunings confirms that this effect is not due to absorption, 
but quantum interference between the intra-cavity field and the ions. b. High resolution transmission spectrum (red curve) 
obtained by scanning a narrow linewidth (~20 KHz) Ti:Sapphire laser over the inhomogeneous line. Green curve is the fit 
using parameters: g— 27rx6 MHz, Th^lOO kHz, Fi n h O m=16.0 GHz, and assuming a Gaussian ion density distribution. The 
green shaded region is the estimated fluctuation in the transmitted laser intensity caused by \/TV statistical variations of the 
ion density. Large fluctuations are expected because the density N is low (few tens), which agrees with the measurement. The 
fluctuations within the inhomogeneous linewidth is noticeably larger than that at far detunings (>25 GHz) and the technical 
background noise (grey area), confirming that they are caused by SFS of the ion spectral density, c. Two traces of the 
transmitted intensities over the same 100 MHz bandwidth near zero detuning at different times. The high degree of correlation 
confirms the static and repeatable nature of SFS. 


probe light power and observing the saturation of the 
ions at increasing intra-cavity photon number (Supple¬ 
mentary S7). The saturation photon number in the nano¬ 
cavity was measured to be (n cav ) = 2 x 10 -5 . 

To better resolve the spectrum, a narrow (^20 kHz) 
Ti:Sapphire laser was scanned across the resonance (See 
Methods) to give the transmitted signals shown in 
Fig. gb, c). A 75 % decrease in transmission was mea¬ 
sured at zero detuning, which corresponds to a collec¬ 
tive cooperativity 77 ^ 1.2. Fitting using a Gaussian 
spectral density distribution (green line) with measured 
parameters g=2nx6 MHz, Fh=100 kHz, ri n h O m=16.0 
GHz, gives a peak ion density of N ~53. Due to the 
statistical fine structure (SFS) of the inhomogeneously 
broadened line, a variation in the transmitted intensity 
(( ST ) 2 ~ ?(5N ) 2 ) due to SN = y/N fluctuations in 
the ion spectral density is expected [31j. This expected 
variation is represented in Fig. [4^b) by the green shaded 
region and shows good agreement with that of the mea¬ 
sured signal. This variation within the inhomogeneously 
broadened line, on which the statistical fluctuations of 
the ion spectral density are imprinted, is significantly 
larger than the spectrometer technical noise (grey area) 
and laser shot noise at far detunings (>25GHz), thus con¬ 
firming that the static SFS in spectral density N(AX) is 
probed. Two traces of the laser scan over the same 100 
MHz bandwidth near zero detuning at different times 
are shown in Fig. [4jc). The high degree of correlation 
reflects the static and repeatable nature of SFS. Notably, 


the current platform would allow detection and control 
of a single ion coupled to the cavity if N <<1 and the 
laser linewidth were narrower than (Supplementary 

S 8 ). 


DISCUSSION 

The results reported in this paper (long optical coher¬ 
ence time, small inhomogeneous broadening, enhanced 
coherent optical control, and resonant probing of cavity- 
coupled REIs) demonstrate REI-based nanophotonics 
as a promising approach for robust and scalable quan¬ 
tum photonic networks integrating memories and sin¬ 
gle REI qubits. Single photon rates exceeding 1 MHz 
can be achieved with single REIs in nano-cavities with 
Q/V ^ 10 4 - 10 5 (V is normalized to (A/n) 3 ), and the in¬ 
homogeneous broadening allows for frequency multiplex¬ 
ing of multiple REIs. To use the interface as an optical 
quantum memory, efficient optical pumping into the long- 
lived Zeeman level needs to be demonstrated. Bulk REI 
quantum memories already boast high storage efficiency 
[32j with multi-mode capacity [3]. Their implementa¬ 
tions in our nanophotonic platform open the possibility 
of multiplexed systems for on-chip quantum repeaters. 
For Nd, high-fidelity storage of entanglement based on 
atomic frequency comb (AFC) has been demonstrated 
[33] [34] . With cavity impedance matching |35], unit stor¬ 
age efficiency is achievable with a mesoscopic ensemble of 
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cavity-coupled ions. Meanwhile, long-lived nuclear spin 
coherence of 9 ms in 145 Nd 36 bodes well for spin-wave 
quantum memories using our nano-resonators. These de¬ 
vices can be further coupled to superconducting or op¬ 
tomechanical devices to enable hybrid quantum systems 
[4]. Furthermore, the technology can be readily trans¬ 
ferred to other wavelengths, such as 1.5 fi m for telecom 
quantum memories using Er 3+ :YSO or 590 nm for long- 
haul quantum hard drives using Eu 3+ :YSO [16]. 

METHODS 

YSO nano-resonator design and fabrication. The 

nano-beam has an equilateral triangular cross-section 
with each side of 780 nm. This geometry allows a cir¬ 
cular fundamental mode field that can be efficiently cou¬ 
pled with a free space laser beam. The cavity is formed 
by forty equally spaced grooves of lattice constant 340 
nm on the nano-beam, except for a defect introduced 
at the center by perturbing the lattice constant. The 
depth of the grooves is 65% of the beam height. The 
triangular nano-beam resonator was fabricated using fo¬ 
cused ion beam (FIB) milling followed by wet etching of 
Ga+ contaminated YSO in diluted (10%) hydrochroric 
acid. An ion beam of 20 kV, 0.2 nA was used to fabri¬ 
cate the suspended nano-beam waveguide by milling at 
30° angle with respect to the crystal surface normal. A 
small ion beam of 23 pA was then used to accurately 
pattern the grooves on top of the nano-beam. Limited 
by the finite width of the FIB beam, the side-walls of the 
grooves in the actual device were not vertical, but had 
an angle of 6°. This leads to a degraded theoretical Q of 
5.0 x 10 4 . We were able to reproducibly fabricate arrays 
of resonators (up to six) in a batch (Supplementary In¬ 
formation SI), with all the devices measuring resonances 
close to 883 nm and quality factors varying from 1,100 
to 10,000. 

Experimental setup for the photon echo measure¬ 
ments. A 500 mT external magnetic field was ap¬ 
plied at 0=135° relative to the crystal Di axis using 
a pair of permanent magnets (see Fig. [ljb)) . The ir/2 
and 7r Gaussian pulses were generated by amplitude- 
modulating the Ti:Sapphire laser with two acousto-optic 
modulators (AOM) in series, with each in a double-pass 
configuration. The two pulse widths were 250 ns and 400 
ns at a repetition rate of 1 kHz. The average (peak) 
power of the excitation pulses was 210 nW (320 /iW) 
measured after the objective lens. The extinction ratio 
between the on and off level of the pulses was ^120 dB, 
ensuring sufficient signal-to-noise ratio for detecting echo 
photons using a Si single-photon counter (Perkin Elmer 
SPCM). A third shutter AOM in single-pass configura¬ 
tion was inserted just before the photon counter to block 
the strong excitation pulses from saturating the detector. 
The extinction ratio of this shutter AOM was 30 dB. 


High-resolution laser spectroscopy on cavity- 
coupled Nd 3+ ions. For the cavity transmission exper¬ 
iments, the Ti:Sapphire laser (M Squared SolsTiS) was 
continuously scanned at a rate of 10 MHz per second. 
The high-sensitivity CCD camera in the spectrometer 
(Princeton Instruments PIXIS) registers the transmitted 
photon counts intermittently at an adjustable frame rate 
(frames per second (fps)) with an exposure time 0.01 s 
for each frame. Therefore one exposure corresponds to a 
spectral width of 10 MHz x 0.01=100 kHz scanned by the 
laser, which is equal to the homogeneous linewidth of the 
0.2% doped sample. The long term drift of the laser is 10 
MHz per hour, so the drift during each exposure should 
be inconsequential. Each data point in Fig. |IJb,c) rep¬ 
resents the photon count collected in one camera expo¬ 
sure, corresponding to the signal contributed by the ions 
within one homogeneous linewidth. The frame rate was 
0.1 fps for the coarser scan in Fig. Qb), corresponding 
to a spectral interval of 100 MHz between two adjacent 
data points. The frame rate was 8.2 fps for the fine scan 
in Fig. [4jc), with a spectral interval ~ 1 MHz. Each 
data point was obtained with one scan. Several scans 
at different spectral regions were performed and stitched 
together to cover the entire bandwidth in Fig. j4^b). 
Estimation of the branching ratio. The measured 
optical depth of a 15-/im-long nano-beam resonator at 
3.8 K was d= 0.1, from which we deduce an oscillator 
strength of / = 6.5 x 10 -7 and a spontaneous emission 
rate of this transition to be 7883 = 1/t 883=1/5.6 ms [24]. 
With a measured bulk medium lifetime To=250 fis , the 
branching ratio was thus estimated to be tq/t 883 ~ 4.5%. 
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SUPPLEMENTARY INFORMATION 


SI. Fabrication and characterization of arrays of YSO nano-beam resonators 


The YSO nano-resonators were fabricated in batch with careful focused ion beam (FIB) alignment and drift com¬ 
pensation. Figure SI (a) shows a SEM image of an array of six Nd:YSO nano-resonators. All of the devices have 
resonance modes near the designed 883 nm wavelength, shown in the transmission spectrum in Fig. |Sl[b). The color 
of each spectrum maps to that of the device in Fig. SI (a). Measured quality factors in this batch range from 1,500 
to 4,400. The spread of resonance wavelengths is about 10 nm, indicating the robustness of this fabrication process. 



b 



875 880 885 

A(nm) 


FIG. SI. An array of nano-beam resonators fabricated in Nd:YSO. a, Scanning electron microscope image of 6 devices, 
b, Measured resonance modes (color matched to the corresponding device) near the 883 nm transition of Nd:YSO (grey line). 


S2. Measurements of inhomogeneous linewidths for cavity-coupled ions 

The inhomogeneous linewidth for the ions coupled to the 0.2% nano-cavity was measured to be Ai n h O m=16.0 GHz 
from the dipole-induced transparency signal in Fig. 4(b). For the 0.003% cavity, ion density was too low for a 
similar measurement of the inhomogeneously broadened distribution. Instead, the linewidth of photoluminescence 
(PL) from the cavity was measured with a high resolution spectrometer, and A in hom=5.9 GHz was estimated by 
deconvolving the PL signal with the minimally resolvable linewidth of the spectrometer. For both doping levels, the 
same inhomogeneous linewidths were measured from the bulk via absorption spectroscopy. The agreement between 
the inhomogeneous linewidth of the cavity-coupled Nd ions and the bulk confirms the excellent spectral stability of 
REIs when embedded in nanophotonic resonators. 


S3. Requirement on the inhomogeneous linewidths for scalable QLMIs 

We consider a network of QLMIs each being a nano-cavity coupled to ensembles of emitters with inhomogeneous 
linewidth Ai n hom- Efficient QLMIs require the emitters to emit photons dominantly into the cavity mode. The 
cavity photons in a single spatial mode could then be efficiently coupled to waveguides or fibres for routing to other 
QLMIs operating at the same frequency. The probability of an emitter to emit a photon into the cavity mode is 
/3F/( 1 + (F — l)/3), where (3 is the branching ratio of the dipole transition, and F is the Purcell factor in Eq. [3j 
Assuming 99% of the dipole emission into the cavity and typical branching ratio of (3 ~ 10%, the required Purcell 
factor should be ^1,000. Considering a photonic crystal nano-cavity with a small mode volume of l(A/n) 3 , F ~ 1, 000 
corresponds to a quality factor Q ~ lx 10 4 and a cavity linewidth k ^30 GHz (for the 883 nm transition). Thus, 
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for scalability, the emitters and cavities need to be aligned within «10 GHz (order of magnitude). This limits the 
inhomogeneous broadening of the dipole ensembles to be <10 GHz for implementing robust and scalable QLMIs. This 
condition is satisfied by most REI transitions. 


S4. Calculation of ensemble averaged Purcell enhancement factor 


The spontaneous emission rate of a dipole coupled to a nano-resonator is enhanced relative to the bulk medium, by 
the factor ~ 1 + (3F pQ, where /3 is the branching ratio of the transition, and F is given by [2], 

i-'o-)-/' v J i 


F F ° av \ |E ma x|H / 1 + 4Q 2 (A/A cav -1) 2 (2) 

where /i is the dipole moment, E(r) is the local electric field at the emitter location r, A cav is the cavity resonant 
wavelength, A is the emitter wavelength, and |E max | is the maximum electric field in the resonator. For a dipole that 
is resonant with the cavity and ideally positioned and oriented with respect to the maximum cavity field, 

3 / A rav \ 3 Q 


(3) 

We consider an ensemble of Nd ions uniformly distributed inside the YSO cavity. The enhancement of the emission 
from the ensemble can be estimated by averaging F (Eq. [2| over the entire population of Nd ions in the cavity. Based 
on the 3 dimensional field profile in Fig. 1(b), the mode volume 1.65(A/n) 3 and Q=4,400, we numerically calculate 
this averaged Purcell factor to be 45 when the cavity is resonant with the transition. 

If the emission rate for uncoupled Nd ions is 1/to = 1/ts83 + l/^other 5 in the coupled case the rate becomes 
l/r c = (1T -E) / t 883 Y1 /t other •> where 1 /ts 83 and l/r 0 ther are the spontaneous emission rates into the 883 nm transition 
and other 4f-4f transitions, respectively. The Purcell factor is then experimentally extracted as F = (to/t c — 1)//?, 
where /3 is the branching ratio of the 883 nm line. Based on the measured branching ratio /3= 4.5%, the observed 
change in lifetimes leads to an ensemble averaged Purcell factor ~ 42, which matches well with the calculated value. 
Furthermore, the averaged value of 42 means the expected Purcell enhancement for an ideally positioned and oriented 
Nd dipole is F cav ~ 200. 


(/) 
H—* 

c 

o 

O 



FIG. S2. Lifetime measurements for a 0.003% doped Nd:YSO resonator. Lifetime changes from 289 /is when the 
cavity of Q=l,500 is detuned by A=0.4 nm, to 200 /is when the cavity is resonant with the ions. 
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S5. Purcell enhancement in 0.003% doped Nd:YSO nano-resonators 

A 0.003% doped Nd:YS0 nano-resonator was fabricated, measuring a resonance mode at 879 nm with quality factor 
Q= 1,500. Spontaneous emission rate enhancement in this cavity was estimated from lifetime measurements in the 
same way as for the 0.2% cavity. A change of lifetime from 289 fis when the cavity is detuned by A =0.4 nm, to 195 
fis at resonance ((Figure [Si]) gives rise to an ensemble averaged Purcell factor 10. Note that a longer T\ = 290 /is 
in the low density Nd:YSO bulk sample (T\ = 300 /is reported in [3]) yields a slightly larger branching ratio /3 ~ 5.4%. 
The longer T\ is most likely due to weaker dipole-dipole interactions in 0.003% doped sample. 


S6. Spectral diffusion of Nd 3+ ions coupled to the nano-cavities 

The dynamic coherence properties of the cavity-coupled Nd ions were investigated by three-pulse photon echoes ((tt/2 
- 7r/2 - 7r/2)) (Figure |S3^a) )that gives information about the spectral diffusion on time scales up to T\ [4] . The 
third pulse, delayed by a time T w after second pulse, is diffracted on the spectral grating from the first two pulses 
and produces an echo. Spectral diffusion - frequency shifts of the optical transition due to the fluctuating rare earth 
environment - gradually erases the grating during T w , and causes faster echo decays thus broadening of the effective 
linewidth T e ff. Linearly increasing T e ff = Th + RT W at a rate of R = 380 Hz /is -1 was measured for the 0.003% 
cavity and 6.1 kHz /is -1 for the 0.2% cavity (Figure [S3^b,c)). Higher spectral diffusion is expected for higher doping 
because of stronger dipole-dipole interaction between Nd ions. Nevertheless, the measured linewidth broadening is 
much smaller than our Rabi frequency (~6 MHz in Fig. 3(c)). This indicates the coupled ions, either singles or 
ensembles, can be optical addressed repeatedly up to 10s of /is, which is desirable for optical quantum information 
processing. 


a 



T 

w 



time 


te/ 2 echo 


b c 




FIG. S3. Spectral diffusion of cavity-coupled ions a, Three-pulse photon echo sequence (tt/2 - 7t/ 2 - 7r/2). b-c, T w 
dependent broadening of the effective linewidths of ions coupled to the 0.2% (b) and 0.003% (c) doped cavities, measured via 
three-pulse photon echoes. Linear fits indicate the spectral diffusion rates R for both doping concentrations. 
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S7. REI-controlled cavity transmission and saturation of the coupled ions 

Figure |S4| plots the on resonance transmission at zero detuning as a function of the average photon number in the 
cavity (n cav ). (n cav ) was estimated from the input probe laser power Pi n (measured after the objective), the coupler 
efficiency 77 , and cavity coupling rate (k/ 2) as (n cav ) = Black curve is the theoretical calculation using the 

Quantum toolbox [5], which shows close agreement with the experiment. 


CD 



FIG. S4. Controlled cavity transmission versus intra-cavity photon number for probe laser at zero detuning. The 

cavity transmission is normalized by empty-cavity transmission at the zero detuning. Saturation occurs at (n cav ) ~ 2 x 10 -5 . 
The inset shows a transmission spectrum at onset of saturation with a normalized transmission at the dip ~54%. 


S8. Towards detection and control of single REI ions coupled to the nano-cavity 

The measurement of 7V(AA) indicate that this system can be used to detect and control a single ion coupled to the 
cavity - a key ingredient for realizing quantum networks interconnecting multiple quantum bits encoded in individual 
REI ions. In the 0.003% low density devices, we estimated a peak ion density of 7V=0.07 per Fh=l/ 7 rX 2 = 3.1 kHz. 
Correspondingly, the single ion cooperativity of 77 =1.6 can be attained with the same cavity Q=4,400, V = 1.65(A/n) 3 , 
and g = 27 tx 10 MHz (typical for REI transitions and we assume the ion is positioned at maximum cavity field). 
Simulation using Quantum Optics Toolbox [5] yields a transmission dip >80% due to a single Nd ion, as shown 
in Fig. |S5| The main technical challenge to detect single ion in this cavity system is the requirement of a highly 
stabilized laser, with linewidth <1 kHz and minimal long term drift, for scanning the single ion spectrum, which 
should be attainable with state of the art laser spectroscopy technology. 


faraon@caltech.edu 
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FIG. S5. Cavity transmission when coupled to a single Nd +3 ion. The simulation uses Quantum Optics Toolbox ^ 
with parameters Q=4,400, Th=3.1 kHz, and g= 27rxl0 MHz. The transmission dip has a full-width at half-maximum (FWHM) 
of ~15 kHz. 
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